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Abstract 

Large  eddy  simulation  (LES)  is  conducted  of  the  Sandia  Flame  D  [Proc.  Combust.  Inst.  27  (1998)  1087, 
Sandia  National  Laboratories  (2004)],  which  is  a  turbulent  piloted  nonpremixed  methane  jet  flame.  The 
subgrid  scale  (SGS)  closure  is  based  on  the  scalar  filtered  mass  density  function  (SFMDF)  methodology 
[J.  Fluid  Mech.  401  (1999)  85].  The  SFMDF  is  basically  the  mass  weighted  probability  density  function 
(PDF)  of  the  SGS  scalar  quantities  [Turbulent  Flows  (2000)].  For  this  flame  (which  exhibits  little  local 
extinction),  a  simple  flamelet  model  is  used  to  relate  the  instantaneous  composition  to  the  mixture  fraction. 
The  modelled  SFMDF  transport  equation  is  solved  by  a  hybrid  finite-difference/Monte  Carlo  scheme.  This 
is  the  first  LES  of  a  realistic  turbulent  flame  using  the  transported  PDF  method  as  the  SGS  closure.  The 
results  via  this  method  capture  important  features  of  the  flame  as  observed  experimentally. 
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1.  Introduction 

There  have  been  significant  recent  develop¬ 
ments  of  subgrid  scale  (SGS)  closures  for  large 
eddy  simulation  (LES)  of  turbulent  reacting  flows. 
Several  recent  reviews  are  available  [4-12].  One 
such  closure  is  via  the  filtered  density  function 
(FDF)  methodology,  first  introduced  by  Pope 
[13].  This  is  the  counterpart  of  the  probability 
density  function  (PDF)  method  that  has  proven 
quite  effective  in  Reynolds  averaged  simulations 
(RAS)  [4,14].  This  success  is  due  to  the  inherent 
property  of  the  PDF  to  provide  complete  statisti- 
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cal  information  about  the  variables.  Due  to  this 
property,  the  FDF  offers  the  ability  to  treat  fi¬ 
nite-rate  chemistry  and  the  turbulence-chemistry 
interactions.  In  comparison  to  RAS,  the  LES/ 
FDF  methodology  provides  a  more  detailed  and 
reliable  prediction  of  turbulent  reacting  flows 
and  is  better  suited  to  account  for  the  large  scale 
unsteady  phenomena,  which  are  prevalent  in  com¬ 
bustion  devices.  The  scalar  FDF  (SFDF)  is  con¬ 
sidered  by  many  investigators  [15-19].  Its 
extension  to  account  for  variable  density  flows  is 
via  the  scalar  filtered  mass  density  function 
(SFMDF)  [3,20].  The  velocity  FDF  (VFDF)  is 
introduced  by  Gicquel  et  al.  [21],  and  the  veloci¬ 
ty-scalar  FDF  (VSFDF)  by  Sheikhi  et  al.  [22].  A 
recent  review  on  the  state-of-progress  on  LES/ 
FDF  is  provided  by  Givi  [12].  The  outcome  of 
these  investigations  has  been  encouraging, 
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warranting  further  extension  and  application  of 
this  methodology  in  turbulent  combustion. 

In  this  work,  we  employ  the  FDF  method  for 
prediction  of  the  piloted  jet  flame  studied  in  the 
experiments  of  the  Combustion  Research  Facility 
at  the  Sandia  National  Laboratories  [1 ,23,24].  This 
flame  has  been  the  subject  of  broad  investigations 
by  other  computational/modelling  methodologies 
[2,25-28].  In  the  experiments,  three  turbulent 
flames  are  considered:  Flames  D,  E,  and  F.  The 
geometrical  configuration  in  these  flames  is  the 
same,  but  the  jet  inlet  velocity  is  varied.  In  Flame 
D,  the  fuel  jet  velocity  is  the  lowest,  and  the  flame 
is  close  to  equilibrium.  The  jet  velocity  increases 
from  flames  D  to  E  to  F,  with  noticeable  non-equi¬ 
librium  effects  in  the  latter  two.  Flame  D  is  consid¬ 
ered  in  this  work.  The  objective  was  to  assess  the 
predictive  capability  of  the  LES/FDF  methodol¬ 
ogy  in  capturing  the  flow  field  and  scalar  mixing. 
This  is  a  necessary  step  before  consideration  of 
the  non-equilibrium  flames  (E  and  F). 


2.  Formulation 

Sandia  Flame  D  consists  of  a  main  jet  with  a 
mixture  of  25%  methane  and  75%  air  by  volume. 
The  nozzle  is  placed  in  a  coflow  of  air,  and  the 
flame  is  stabilized  by  a  substantial  pilot.  The  Rey¬ 
nolds  number  for  the  main  jet  is  Re  =  22400  based 
on  the  nozzle  diameter  D  =  1.2  mm  and  the  bulk 
jet  velocity  49.6  m/s.  All  of  the  details  regarding 
this  flame  are  provided  and  regularly  updated  in 
the  web  site  [2].  The  methane-air  reaction  mecha¬ 
nism,  as  occurs  in  this  flame,  is  taken  into  account 
via  the  “flamelet”  model.  This  model  considers  a 
laminar,  one-dimensional  counterflow  (opposed 
jet)  flame  configuration  [8].  The  detailed  kinetics 
mechanism  of  the  Gas  Research  Institute 
(GRI2.11)  [29]  is  employed  to  describe  combus¬ 
tion.  At  low  strain  rates,  the  flame  is  close  to 
equilibrium.  Thus,  all  of  the  thermo -chemical 
variables  are  determined  by  the  “mixture 
fraction.” 

Formulation  is  based  on  the  compressible  form 
of  the  continuity,  Navier-Stokes,  energy  (enthal¬ 
py)  and  mixture  fraction  conservation  equations 
in  a  low  Mach  number  flow  [30].  These  equations 
govern  the  space  (x  =  xh  i  =  1,2,3)  and  time  {t) 
variations  of  the  fluid  density  p(x,  t),  the  velocity 
vector  u  =  wf-(x,  /),  the  pressure  p(x,  t ),  the  specific 
enthalpy  h(x,t ),  and  the  mixture  fraction  £(x,  t). 
We  employ  Fourier’s  law  of  heat  conduction, 
Fick’s  law  of  diffusion,  and  we  assume  unity  Lewis 
number.  The  molecular  viscosity  /a  increases  with 
temperature  (7)  to  the  power  of  0.7.  The  magni¬ 
tude  of  the  molecular  Schmidt  (and  Prandtl)  num¬ 
ber  is  Sc  =  0.75. 

For  LES,  we  consider  a  spatially  and  tempo¬ 
rally  invariant  positive  filter  Q(x!  -  x)  of  width 
Ag.  The  filtered  value  of  the  variable  Q{x,  t)  is 


denoted  by  (Q(x,t))e,  and  its  mass-weighted  fil¬ 
tered  value  is  denoted  by  (Q(x,  t))L  =  (pQ)e/(p)e- 
The  SFMDF  of  the  mixture  fraction  is  denoted 
by  Tl{ i/qx,  f),  where  0  denotes  the  (probability) 
composition  domain  of  the  mixture  fraction. 
The  SFMDF  accounts  for  SGS  statistics  of  only 
the  scalar  field.  Thus,  the  closure  of  the  SGS 
hydrodynamics  must  be  provided  by  other 
means.  Here,  we  employ  a  well-established  gra¬ 
dient  diffusion  model,  in  which  the  SGS  dy¬ 
namic  viscosity  vt  is  modelled  by  the  MKEV 
(modified  kinetic  energy  viscosity)  [3,31]  model. 
In  the  transport  equation  governing  the 
SFMDF,  the  effects  of  SGS  convection  are  also 
modelled  by  an  analogous  gradient  diffusion 
model.  The  SGS  mixing  is  closed  via  the  least 
mean  square  estimation  (LMSE)  model  [32,33]. 
These  are  described  in  detail  in  our  previous  pa¬ 
pers  on  SFDF  [3,17].  The  final,  modelled 
SFMDF  transport  equation  reads 

+  mW-®L)FL\.  (1) 

Here,  y  =  /ll/ Sc  is  the  molecular  diffusivity, 
yt  =  is  the  SGS  diffusivity,  and  Sct  =  0.75 
is  the  SGS  Schmidt  number.  With  the  closure 
as  such,  the  modelled  scalar  flux  is  consistent 
with  those  in  most  other  (non-FDF)  methods 
[34,35] 

(2) 

The  term  Qm  is  the  SGS  mixing  frequency  (as  re¬ 
quired  in  LMSE)  and  is  modelled  by  [3,17] 
4n  =  Cn(y  +  7,)/((p)e^2g),  with  CQ  =  8. 

The  relation  between  the  thermo-chemical 
variables  (denoted  by  the  array  (/){x,t ))  and  the 
mixture  fraction  is  provided  by  the  flamelet  li¬ 
brary.  In  the  context  of  the  opposed  jet  flame, 
the  library  provides:  0  =  0(£,a)  where  a  denotes 
the  strain  rate.  For  a  fixed  value  of  a ,  the  SGS  sta¬ 
tistics  of  the  thermo-chemical  variables  are  deter¬ 
mined  from  the  SFMDF 

<0(X,  t))L  =  A  J  (/>(£)  fL(Z,  X,  t)  d£.  (3) 

A  hybrid  finite-difference  (FD)/Monte  Carlo 
(MC)  scheme  is  employed  to  solve  the  coupled 
set  of  the  SFMDF  equation  and  the  filtered 
hydrodynamic  equations.  In  this  scheme,  the  do¬ 
main  is  discretized  by  FD  grid  points,  and  the 
SFMDF  is  represented  via  an  ensemble  of  MC 
particles  [36].  All  of  the  hydrodynamic  variables 
are  determined  on  the  FD  points.  A  fourth  order 
compact  scheme  is  used  for  FD  discretization  of 
the  compressible  flow  equations,  as  described  in 
[37,38].  Transport  of  the  MC  particles  and  the 
change  in  their  properties  are  described  by  a  set 
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of  stochastic  differential  equations  (SDEs)  similar 
to  those  describing  a  diffusion  process  [39].  The 
MC  particles  undergo  motion  in  physical  space 
by  convection  due  to  the  filtered  flow  velocity 
and  diffusion  due  to  molecular  and  SGS  diffusivi- 
ties.  The  compositional  value  of  each  particle  is 
changed  due  to  SGS  mixing.  These  are  described 
by  the  SDEs: 


(P)e  fai  \ 

d  t 

\J2(y  +  y,)/(p)tdwi-. 

(4) 

om[e 

(5) 

where  xf  and  £+  denote  Lagrangian  position  and 
the  composition  (mixture  fraction),  respectively. 
The  term  w  denotes  the  Wiener  process  [40].  The 
Fokker-Planck  equation  corresponding  to  Eqs. 
(4),  (5)  is  equivalent  to  Eq.  (1).  Thus,  the  solution 
of  these  SDEs  represents  the  FMDF  in  the  prob¬ 
abilistic  sense. 

To  understand  the  operational  procedure,  the 
elements  of  the  computation  are  shown  in  Fig.  1. 
For  simplicity,  a  two-dimensional  domain  is 
shown  with  fixed  grid  spacing  of  size  of  A.  The 
MC  particles  are  distributed  randomly  and  are 
free  to  move  within  the  domain  as  governed  by 
Eq.  (4).  This  transport  is  Lagrangian,  thus  the 
solution  is  free  of  constraints  associated  with 
typical  simulation  of  convection  on  fixed  grid 
points.  Statistical  information,  e.g.,  filtered  val¬ 
ues,  at  any  point  is  obtained  by  considering  an 
ensemble  of  NE  computational  particles  residing 


Fig.  1.  Elements  of  computation  as  used  in  a  typical 
LES/FDF.  Solid  squares  denote  the  FD  points,  and  the 
circles  denote  the  MC  particles.  Also  shown  are  three 
different  ensemble  domains. 


within  an  ensemble  domain  of  side  length  AE 
centered  around  the  points.  For  reliable  statistics 
with  minimal  numerical  dispersion,  it  is  desired 
to  minimize  the  size  of  ensemble  domain  and 
maximize  the  number  of  MC  particles  [14].  In 
this  way,  the  ensemble  statistics  would  tend  to 
the  desired  filtered  values.  Transfer  of  informa¬ 
tion  from  the  grid  points  to  the  MC  particles  is 
accomplished  via  interpolation.  The  transfer  of 
information  from  the  particles  to  the  grid  points 
is  accomplished  via  ensemble  averaging.  In  the 
hybrid  scheme,  some  of  the  filtered  quantities 
are  obtained  by  MC,  some  by  FD,  and  some 
by  both.  That  is,  there  is  a  “redundancy”  in 
determination  of  some  of  the  quantities.  This 
redundancy  is  very  useful  in  monitoring  the  con¬ 
sistency  of  the  simulated  results  [41].  In  general, 
all  of  the  equations  for  the  filtered  quantities 
can  be  solved  by  FD,  where  all  of  the  unclosed 
terms  are  evaluated  by  MC.  This  process  can 
be  done  at  any  filtered  moment  level  [21].  Here, 
the  filtered  values  of  the  mixture  fraction  and 
the  temperature  are  used  to  establish  consistency. 
Some  results  in  this  regard  are  presented  in  the 
next  section.  For  more  details  we  refer  to 
[3,17,21,22].  In  addition,  the  computational  accu¬ 
racy  of  the  methodology  is  established  by  exam¬ 
ining  both  the  statistical  and  the  dispersion 
errors  [42].  In  doing  so,  the  correlation  of  the 
fluid  density  with  distribution  of  the  MC  parti¬ 
cles,  the  size  of  the  ensemble  domain,  the  number 
of  particles  within  this  domain,  and  the  global 
distribution  of  the  particles  are  monitored  in  a 
manner  similar  to  those  reported  in  our  previous 
work  [3,17,21,22]. 


3.  Results 

The  flow  variables  at  the  inflow  are  set  the 
same  as  those  in  the  experiments,  including  the  in¬ 
let  profiles  of  the  velocity  and  the  mixture  frac¬ 
tion.  The  flow  is  excited  by  superimposing 
oscillating  axisymmetric  and  helical  perturbations 
at  the  inflow.  The  procedure  is  similar  to  that  of 
Danaila  and  Boersma  [43],  but  the  amplitude  of 
forcing  is  set  in  such  a  way  as  to  match  the  exper¬ 
imentally  measured  turbulent  intensity  of  the 
streamwise  velocity  at  the  inlet.  Simulations  are 
conducted  on  a  three-dimensional  Cartesian  mesh 
with  uniform  spacings  in  each  of  the  three  direc¬ 
tions.  Standard  characteristic  boundary  condi¬ 
tions  [44]  are  implemented  in  all  of  the  FD 
simulations.  The  computational  domain  spans  a 
region  of  18D  x  lOD  x  lOD  in  streamwise  (x), 
and  the  two  lateral  (y,z)  directions,  respectively. 
The  number  of  grid  points  is  91x101x101  in 
the  x,  y,  and  z  directions,  respectively.  The  filter 
size  is  set  equal  to  Ag  =  2(AxAyAz)(1/3)  where 
Ax,  Ay,  and  A z  denote  the  grid  spacings  in  the  cor¬ 
responding  directions.  The  size  of  the  ensemble 
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domain  for  evaluations  of  the  filtered  values  from 
the  MC  solver  is  equal  to  the  filter  size.  The  par¬ 
ticles  are  supplied  in  the  inlet  region  and  are  free 
to  move  within  the  domain  due  to  combined  ac¬ 
tions  of  convection  and  diffusion  (molecular  and 
SGS).  There  are  approximately  48  MC  particles 
at  each  grid  point.  As  per  results  of  extensive  pre¬ 
vious  studies  [3,17,21,22],  this  is  sufficient  to  yield 
an  excellent  statistical  accuracy  with  minimal  dis¬ 
persion  errors.  In  total,  there  are  about  3.4  million 
MC  particles  within  the  domain  at  all  times.  The 
simulation  results  are  monitored  to  ensure  that 
the  particles  fully  encompass  and  extend  well  be¬ 
yond  regions  of  non-zero  vorticity  and  reaction. 

The  flamelet  table  at  a  strain  rate  of  a  =  100 1  /s 
is  used  to  relate  the  thermo -chemical  variables  to 
the  mixture  fraction.  This  value  is  consistent  with 
those  used  in  previous  PDF  predictions  of  Sandia 
flames  [45].  First,  the  consistency  and  accuracy  of 
the  computations  are  assessed.  Next,  the  overall 
predictive  capability  of  LES/SFMDF  is  demon¬ 
strated  by  comparing  the  flow  statistics  with  the 
Sandia  data.  These  statistics  are  obtained  by 
long-time  averaging  of  the  filtered  field  during  6 
flowthrough  times.  A  total  of  50,000  samples  of 
each  of  the  variables  are  collected  in  this  recording 
period.  The  notations  Q  and  RMS(g)  denote, 
respectively,  the  time-averaged  mean  and  root 


mean  square  values  of  the  variable  Q.  The  sim¬ 
plest  consistency  check  is  via  flow  visualization. 
For  example,  Fig.  2  shows  the  instantaneous  con¬ 
tours  of  the  filtered  mixture  fraction  field  as  ob¬ 
tained  by  the  FD  and  the  MC  methods.  The 
central  jet  lies  in  the  middle  along  the  axial  coor¬ 
dinate,  surrounded  by  a  pilot  where  the  tempera¬ 
ture  is  the  highest  and  encircled  by  the  air  coflow. 
The  region  close  to  inlet  is  dominated  by  the 
molecular  diffusion,  and  the  jet  exhibits  a  lami- 
nar-like  behavior.  Farther  downstream,  the 
growth  of  perturbations  is  manifested  by  the  for¬ 
mation  of  large-scale  coherent  vortices.  The  up¬ 
stream  feedback  from  the  vortices  created 
initially  triggers  further  self-sustaining  vortex  rol¬ 
lup,  and  subsequent  pairing  and  coalescence  of 
neighboring  vortices  [46,47].  Due  to  the  presence 
of  helical  instabilities,  the  instantaneous  flow  is 
asymmetric.  The  similarity  of  the  results  in  the 
two  figures  is  observed  at  all  other  times  and  is 
also  observed  for  the  temperature  (not  shown). 
This  consistency  is  further  assessed  via  establish¬ 
ing  identical  statistics  of  the  redundant  quantities 
as  generated  by  the  two  methods.  The  capability 
of  the  method  in  predicting  the  hydrodynamics 
field  is  demonstrated  by  examining  some  of  the 
(reported)  flow  statistics.  The  centerline  mean 
and  RMS  values  of  the  axial  velocity  are  com- 


Fig.  2.  Instantaneous  contours  of  the  filtered  values  of  the  mixture  fraction  as  obtained  by  FD  and  MC  simulations. 
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pared  to  experimental  data  in  Fig.  3.  This  figure 
indicates  that  the  flow  is  adequately  excited,  and 
the  predicted  results  are  in  good  agreement  with 
data. 

The  statistics  of  the  thermo-chemical  variables 
are  also  compared  with  corresponding  data.  The 
measurements  are  exhaustive;  here,  only  some 
sample  results  are  shown.  The  radial 
(r  =  \Jz2  +  y2)  distribution  of  the  mixture  fraction 
is  shown  to  compare  well  with  data  (Fig.  4A). 
Similar  agreement  is  observed  at  other  streamwise 
locations.  The  mean  temperature  values  in  Fig.  4B 
indicate  over-prediction  on  the  rich  side.  This  is 
due  to  premixing  of  methane  with  air  as  also  indi¬ 
cated  previously  [27].  The  “resolved”  RMS  values 
of  the  mixture  fraction  and  temperature  are  in 
good  agreement  with  data.  However,  the  “total” 
RMS  values,  including  the  contributions  of  both 
the  resolved  and  the  SGS  fields,  are  higher  than 
values  reported  experimentally.  The  contribution 
of  the  SGS  to  the  total  scalar  energy  is  about 
30%,  which  is  expected  in  LES.  The  higher  values 
of  the  total  RMS  values  as  predicted  by  LES/FDF 
are  not  due  to  MC  numerical  dispersion  because 
the  FD  results  do  indeed  yield  the  same  values. 
The  level  of  SGS  variance  can  be  decreased  by 
increasing  the  magnitude  of  CQ.  However,  this 
would  not  alter  the  total  RMS  values  significantly. 
It  is  possible  that  some  contribution  of  this  vari¬ 
ance  is  not  included  in  the  measurements  due  to 
finite  probe  size.  Higher  resolution  measurements 
would  determine  the  allocations  of  scalar  variance 
to  the  resolved  and  SGS  fields. 

The  statistics  of  the  mass  fractions  (denoted  by 
Y)  of  several  of  the  species  are  compared  with 
data  in  Fig.  5.  Similar  to  the  temperature  results 
on  the  rich  side,  the  reactants’  mass  fractions  are 
under-predicted,  while  the  products  are  over-pre¬ 
dicted.  The  mean  values  of  the  mass  fractions  of 
the  major  and  the  minor  species  compare  well 
with  experimental  data.  All  of  the  results  indicate 


x/D 

Fig.  3.  The  mean  and  RMS  values  of  the  axial  velocity 
(m/s)  at  the  centerline.  Uc l  denotes  the  mean  axial 
velocity  at  the  centerline,  the  symbols  denote  experi¬ 
mental  data,  and  the  lines  denote  the  predictions. 
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Fig.  4.  Radial  distributions  of  the  mean  and  RMS 
values  at  x/D=  15.  The  symbols  denote  experimental 
data.  The  thick  lines  denote  the  mean  values,  the  thin 
solid  lines  denote  the  resolved  RMS  values  and  the  thin 
dashed  lines  denote  total  RMS  values.  (A)  Mixture 
fraction,  (B)  temperature  (K). 


the  adequacy  of  the  flamelet  table  in  relating  the 
thermo-chemical  variables  to  the  mixture  fraction, 
and  also  the  good  predictive  capability  of  the 
SFMDF  for  this  flame.  The  level  of  agreement 
of  the  RMS  values  of  the  mass  fraction  is  the  same 
as  that  presented  in  Fig.  4.  Finally  in  Fig.  6  the 
PDFs  of  the  resolved  mixture  fraction  as  pre¬ 
dicted  by  LES  are  compared  with  those  measured 
experimentally.  It  is  encouraging  to  observe  excel¬ 
lent  agreements  of  both  the  peak  and  the  range  of 
the  PDF.  Similar  agreements  are  observed  at  all 
the  other  locations. 

The  computational  costs  associated  with  LES/ 
FDF  depends,  obviously,  on  the  parameters  of  the 
simulations.  For  the  case  reported  here,  the  simu¬ 
lations  required  about  llOh  of  CPU  time  on  a 
SUN  Fire  4800  with  6  processors.  This  includes 
the  times  required  for  consistency  tests  and 
ensemble  averaging  of  data.  The  computational 
time  for  LES  without  including  SGS  effects  [17] 
is  about  10-12  times  less.  However,  such  simula¬ 
tions  yield  erroneous  predictions  and  in  many 
cases  lead  to  numerical  instabilities.  For  further 
comparative  assessment  of  the  computational 
requirements  of  the  FDF  in  comparison  to  non- 
FDF  methods,  we  refer  to  previous  work 
[3,17,21,22]. 
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Fig.  5.  Radial  distribution  of  the  mean  values  of  the 
mass  fractions  at  x/D  =15.  The  symbols  denote  exper¬ 
imental  data.  (A)  Solid  line,  CH4;  dashed  line,  02;  (B) 
solid  line,  C02;  dashed  line,  CO. 

4.  Summary  and  concluding  remarks 

The  filtered  density  function  (FDF)  methodol¬ 
ogy  [12,13]  is  now  at  a  stage  that  it  can  be  used 
for  prediction  of  complex  turbulent  reacting  flows. 
This  is  demonstrated  in  this  work  by  utilizing  the 
simplest  form  of  the  FDF  for  large  eddy  simula¬ 
tion  (LES)  of  a  piloted,  nonpremixed,  turbulent, 
and  methane  jet  flame  (Sandia  Flame  D).  For  this 
near-equilibrium  flame,  the  thermo -chemical  vari¬ 
ables  are  related  to  the  mixture  fraction.  This  is 
done  by  construction  of  a  flamelet  library  (in  a 
counter-flow  jet  flame)  in  which  the  chemical  reac¬ 
tion  is  modelled  by  detailed  kinetics  [29].  It  is  use¬ 
ful  to  note  that  the  approach  here  is  fundamentally 
different  from  those  followed  in  previous  flamelet 
based  SGS  models.  In  most  previous  contributions 
[27,48-51],  the  FDF  of  the  mixture  fraction  is  as¬ 
sumed  (e.g.,  /3  or  other  distributions).  Here,  a  mod¬ 
elled  transport  equation  for  the  mass  weighted 
FDF  of  the  mixture  fraction  [3]  is  considered.  This 
equation  is  solved  by  a  hybrid  finite-difference/ 
Monte  Carlo  method.  After  establishing  the  con¬ 
sistency  and  accuracy  of  the  hybrid  solver,  the  pre¬ 
dictive  capability  of  the  overall  scheme  is  assessed 
by  comparison  with  experimental  data.  For  these 
comparisons,  the  ensemble  (long  time  averaged) 
values  of  the  thermo-chemical  variables  are  con¬ 
sidered.  It  is  shown  that  all  of  the  mean  quantities 


<  €  >L 


Fig.  6.  PDF  of  filtered  mixture  fraction  at  x/D  =  15  and 
different  radial  locations.  The  symbols  and  the  thick 
lines  denote  experimental  data  and  LES  predictions, 
respectively.  (A)  r/D  =  0,  (B)  r/D  =  1.67. 

are,  generally,  predicted  well.  The  resolved  RMS 
of  the  variables  also  compare  well  with  data.  How¬ 
ever,  when  the  contribution  of  the  subgrid  scale 
(SGS)  quantities  is  included,  the  experimental  data 
are  over-predicted. 

There  are  two  ways  by  which  this  work  can  be 
continued.  First  is  extension  of  LES/SFMDF  for 
prediction  of  flames  that  experience  extinction 
(such  as  Sandia  Flames  E  and  F)  and/or  re-igni¬ 
tion.  This  would  provide  a  more  definitive  assess¬ 
ment  of  the  predictive  capabilities  of  the  FDF. 
Such  simulations  require  consideration  of  finite- 
rate  chemistry,  as  demonstrated  in  RAS/PDF 
simulations  of  Sandia  flames  [25,26].  Presently,  it 
is  not  computationally  feasible  to  implement  de¬ 
tailed  kinetics  in  such  simulations.  But  implemen¬ 
tation  of  reduced  kinetics  schemes  using  in  situ 
adaptive  tabulation,  such  as  those  used  in  RAS/ 
PDF  [25,26],  is  within  reach.  Second,  it  is  desir¬ 
able  to  implement  the  LES/VSFDF  [22]  for 
prediction  of  this  (or  other  complex)  flame(s). 
In  VSFDF,  the  SGS  convection  appears  in 
a  closed  form,  thus  the  assumption  of  gradient 
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diffusion  can  be  relaxed.  This  is  also  feasible  with¬ 
in  the  near  future.  Accomplishments  of  both  of 
these  tasks  can  be  expedited  by  further  reduction 
of  the  computational  requirements  of  the  FDF. 
The  present  work  establishes  the  capability  of 
LES/FDF  for  accurate  prediction  of  complex 
flames,  warranting  its  further  use  for  modelling 
of  even  more  complex  turbulent  reacting  flows. 
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Johannes  Janicka,  Technical  University  Darmstadt,  Ger¬ 
many.  (1)  Could  you  comment  on  the  applied  simulation 
conditions?  Of  interest  are  the  number  of  points,  number 
of  particles,  and  SGS-models  for  scalar  and  momentum. 
(2)  Can  you  comment  on  the  sharp  of  the  SGS-pdf  ? 

Reply.  (1)  The  first  part  of  the  question  is  answered  in 
the  paper. 

(2)  The  structure  of  the  SGS-PDF  depends  on  loca¬ 
tion/time.  No  “general”  comments  can  be  made  about 
its  shape. 


